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ABSTRACT 



The dynamic equations for a Liquid Metal Fast Breeder 
Nuclear Reactor (LMFBR) Power Plant are derived. These non- 
linear differential equations with typical design values for 
the system parameters are solved on a COMCOR Ci 5000 analog 
computer for the characteristic transients of the basic plant 
and the results compared with a recently proposed LMFBR design. 
A digital computer automatic control program is developed 
which uses, in addition to the ideas of derivative and pro- 
portional control, precomputed control trajectories. The 
digital controller design is implemented by applying FORTRAN 
IV programming to the XDS 9300 digital computer. The effec- 
tiveness of the strategy for the automatic controller is 
demonstrated by using various step changes in the load 
demand . 

A special software routine was written to automatically 
scale the analog computer and proved to be quite useful in 
improving the operating efficiency of the hybrid computer 
simulation . 
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I. 



INTRODUCTION 



In view of the enormous increase of world power consumption 
and the fast decrease of the resources which may be used to 
satisfy this growing power demand, it became necessary to 
develop and build nuclear power plants. The reactors used to 
generated the power for these plants are up to now almost 
exclusively thermal reactors using uranium as fissionable 
material. There are 3 fissionable isotopes: U (235), Pu (239), 

and U (233) . U (235) is found in natural uranium and can be 
used in this state or in an enriched form, i.e., its percentage 
is increased beyond that in natural uranium. Pu (239) is 
manufactured by nuclear reaction from U (238) . U (233) is 
similarity manufactored from Th (232) , an isotope which is also 
found in natural uranium. 

The worlds supply of nuclear fuel consists of 25x10^ tons 

g 

of natural uranium, mostly U (238) , and 1x10 tons of thorium, 
mostly Th (232) [10], The natural uranium consists of 99.29% 

U (238) and 0.71% U (235) . Since U (235) is the only naturally 
fissionable isotope the use of nuclear fuels would be limited 

g 

to 0.71% of 25x10 tons. This means that, using only the 
present thermal type reactors with their inefficient nuclear 
economy a shortage of nuclear fuel is lightly to occur by the 
year 2000. 

An increased nuclear economy in an uranium reactor can be 
achieved by using fuel breeding. In this process the nonfis°ion 
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capture of neutrons by a fertile U (238) nucleus results in 
the generation of fissionable Pu (239) . Fertile Th (232) can 
also be converted to fissionable U (233) . If more fissionable 
nuclei are produced than are consumed the reactor in which this 
process takes place is called a Breeder Reactor. If the fission 
process depends on fast neutrons as in the above described 
processes the reactor is a Fast Breeder Reactor (FBR) . 

The control Problem of a FBR is particularly acute because 
of two important facts . First of all due to the high energies 
required of the neutrons in order to produce fission, the power 
per unit volume is much larger than in the presently used 
thermal type reactors fTO] . The high power density results in 
steep temperature gradients, large coolant temperature rises 
and therefore coolant distribution problems. 

The second important fact is the lack of inherent load 
following capability of this type of plant [Xl • This is due 
to the fact that the load fluctuations do not feedback as 
efficiently through the coolant medium as they do in thermal 
type reactors . For these reasons constant automatic super- 
vision and control of a fast breeder reactor is necessary. 

To design a controller for a nuclear power plant it is 
necessary first to obtain a mathematical model. In this report, 

I 

the dynamic equation of a LMFBR are derived. The next step in 
the design procedure was to determine the dynamic behavior 
of the complete plant following selected changes in both the 
controlling inputs (reactivity, valve opening) and in the load. 
Because of the nonlinearity of the system, model simulation 
was required in order to obtain a solution for the dynamic 
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V 




response. From this simulation conclusions were drawn about 
the characteristics which the controller must have and about 
the required control strategy. 

The real-time solution of the dynamic equations is performed 
by the analog computer (COMCOR Ci 5000) in a hybrid linkage 
with the digital computer (XDS 9300) . 

The controller design developed in this report uses the 
digital computer as the controlling element. The economic 
advantage of using a digital controller results from the pros- 
pect of operating a general purpose type computer in the time- 
share mode to perform both the control decisions as well as the 
monitoring data processing and such other related functions 
that are presently being performed by such digital computers. 

A digital computer automatic control program is then 
developed on the basis of the results of the simulated plane 
response. The design strategy includes in addition to the 
ideas of derivative and proportional control, precomputed 
control trajectories. The digital controller design was 
implemented by applying FORTRAN IV programming to the XDS 9300 
digital computer. The effectiveness of the strategy for the 
automatic controller is demonstrated by using various step 
changes in the load demand. 

A special software routine was written to automatically 
scale the analog computer and proved to be quite useful in 
improving the operating efficiency of the hybrid computer 
simulation. 
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II. PLANT DESCRIPTION 



The plant considered in this investigation consists of : 

1) Fast Breeder Nuclear reactor with control rods 

2) Primary sodium coolant loop 

3) Heat exchanger (counterflow) 

4) Secondary sodium coolant loop bypass and bypass 
control valve 

5) Steam generator 

6) Steam loop with turbine and condensor 
Additional components which are found in nuclear reactor 

.power plants such as pumps, intermediate heat exchangers, 
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superheater, and intermediate turbine stages, were not included 
since they would not. have served the purpose of this investi- 
gation. By adding these parts, the basic structure of the 
system would not be changed and would still be valid. 



A flow diagram of the plant is given below. 




Fig. 1. Plant flow diagram. 
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R = Fast Breeder nuclear reactor 
CR = Control rods 
HX = Heat exchanger 
SGen = Steam generator 
T = Turbine 
C = Condenser 
Nal -• Primary sodium loop 
Na2 = Secondary sodium loop 
^0 = Steam circuit 

B = Sodium bypass valve 
TF - Average fuel temperature 
TAV = Average coolant temperature in reactor core 
TCI = Reactor inlet temperature 
THl = Reactor outlet temperature 
TBIl = Inlet temperature heat exchanger primary side 
TB01 - Outlet temperature heat exchanger primary side 
TB1 = Average temperature heat exchanger primary side 
TH2 = Outlet temperature heat exchanger secondary side 
TC2 = Inlet temperature heat exchanger secondary side 
TB2 = Average temperature heat exchanger secondary side 
TBI2 = Inlet temperature steam generator primary side 
TB02 = Outlet temperature steam generator primary side 
TB3 = Average temperature steam generator primary side 
TS = Average temperature steam generator secondary 
side 
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III. MATHEMATICAL MODEL - PLANT DYNAMICS 



A. REACTOR CORE 

1 . Neutron Kinetics 

The neutron kinetic equation with one group approxima 
tion [1] is given as 



dn(t) 

dt 


= n(t) + A c(t) (III.l) 


dc (t ) 
dt 


= | n ( t) - A c ( t ) (III. 2) 


where n(t) = 


neutron density 


6K(t) = 


reactivity 


3 = 


fraction of delayed neutrons 


l = 


neutron lifetime 


A = 


average decay constant 


c ( t) = 


average concentration of delayed neutrons 



The steady state conditions are given by 



<$k(t) 


= 0 


in' = 

e o 


: Xc o 


where n = 

o 


■ steady state neutron density 


c o 


• average steady state concentration of 
delayed neutrons 



13 



When setting 



6K(t) 

a 



K(t) 



= reactivity in dollar units 



the 



n(t) 

n o 

c (t) 
c 

o 

following 



= N(t) 

= D (t ) 

equations are obtained 



= | (K(t) NCt) - N(t) + D (t ) (III. 3) 

= X (N(t) - D(t)) (III. 4) 

The numerical values for a fast reactor for 3 and £ are approx- 
imately given as 

3 s 3 • 10 -3 
£ s 6 • 10~ 7 

Therefore | .5 • 10 4 

If the differential equation for N(t) is brought into the form 
0 = K ( t ) N(t) - N(t) + D (t ) - | 



The last term, being small compared to the other three terms 
may be neglected. Using this approximation the equation for 
the neutron kinetics in the reactor core are then given as 



dD(t) 

dt 

N't) 



A (N(t) 



D ( t ) 
l-K(t) 



D (t ) ) 



(111. 5) 

(111. 6) 
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with the steady state, initial conditions 

N(0) = D (0 ) = 1 (III. 7) 

2 . Reactivity Coefficients 

In addition to the controlled reactivity input applied 
by positioning of the control rods, a nuclear reactor has 
inherent reactivity feedback, produced by thermohydraulic 
effects. The total reactivity is given as 

K (t) = K r (t ) + K p (t) (HI. 8) 

where K_,(t) = controlled reactivity input 
K 

K R (t) = inherent reactivity feedback coefficients 

In a fast reactor the following reactivity coefficients are 
of importance 

a) Doppler coefficient 

b) Fuel axial expansion coefficient 

c) Fuel radial expansion coefficient 

d) Clad temperature coefficient 

e) Sodium temperature coefficient 

a block diagram for the stability model is given in Fig. 2. 
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T 



F 



TAV 



T-CLAD 




Fig . 2 . Block Diagram for reactivity feed back paths . 



The numerical values of the reactivity coefficients are given 



as [2] 




a) 


Doppler coefficient 


b) 


Fuel axial expansions 




coefficient 


c) 


Clad temperature 




coefficient 



- 1.39 cents/ % power 
.32 cents/ % power 

+ .0138 cents/ % power 
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d) Sodium temperature + .0057 cents/ % power 

coefficient 

e) Radial expansion - .0043 cents/ % power 

coefficient 

From the above data it is seen that coefficients c, d, and e 
are several orders of magnitude smaller than a and b. There- 
fore, these reactivity coefficients will be neglected and it 
will be shown in section IV that this does not lead to any 
noticeable inaccuracy. The reactivity feedback coefficient 
may now be expressed as 



k f = k d + k exp 



(III. 9) 



where = doppler coefficient 

K_..„ = fuel axial exoansion coefficient, 
and the equation for the total reactivity is given as 



K T K R + K D + K EXP 



(III. 10) 



If the numerical values for the reactivity coefficients are 
introduced along with the dollar units of reactivity (see 
Appendix A) , the final equation for the reactivity feedback 
is given as 

K m = K - .0014 TF (III. 11) 

1 K. 



with the control block diagram given in Fig. 3. 



17 




Fig. 3. Reactivity Feedback Model 
(Simplified) . 



3 . Heat Generation 

The rate of heat produced per unit volume of fuel is 
called the volumetric thermal source strength and is given 
by (4] 



Q N ■= G • n(t) • a f • cf> (III. 12) 

where G = energy per fission, Mev 

3 

n(t) = number of fissionable fuel nuclei/cm 

Of = Microscopic fission cross section of the fuel 

2 

<j) = neutron flux per cm sec. 

If n is taken as a constant for a certain fuel element, it is 
seen from the above relation, that the generated heat is 
proportional to the product of the fission cross sect:.": 1 . - 
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and the neutron flux <J> . The neutron flux is dependent upon 
the fissionable species used and the energy of fission neutrons, 
which in turn, is dependent upon the reactor temperature. If 
a ^ if determined for a certain reactor, it is seen that the 
heat generated at any point is directly proportional to the 
neutron flux at that point, which is directly proportional to 
the earlier used N (t) . The generated heat may therefore be 
written as 

Q n = A NCt) (III. 13) 

where A = proportionality factor, which is constant if 

Of is assumed to be a constant 

B. HEAT TRANSFER 
1 . Reactor 

a. Fuel to Coolant 

The core of a fast nuclear reactor consists of a 
great number of fuel elements. In a proposed 1 MW Fast 
Breeder Reactor [2] for example, there are 470 fueld rods, each 
fuel section (core length) is 2 feet long and the rod has an 
overall diameter of .25 inches. The fuel is enclosed in a 
cladding of wall thickness of .015 inches. 

Due to the size of the fuel element the radial 
flux distribution within the core may be neglected. The axial 
flux distribution is shown in Figure 4. 
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Coolant 

Fig. 4. Reactore Core and Axial 
Flux Distribution. 



If the axial variation in flux is taken into account, it is 
seen from figure 5 that the generated heat and therefore the 
temperature in the fuel element, cladding and the coolant 
varies in the axial direction of the fuel element. 




Fig. 5. Temperature variation along fuel 
element (center, cladding, 
coolant) . 
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In this investigation, however, since the reactor 
is only a part of the total power plant, these special effects 
will be neglected. Only average temperatures will be 
considered . 




Fig. 6. Temperature profile in 
fuel element. 



Figure 6 shows a typical temperature profile within 
the fuel element, cladding and coolant together with the average 
temperature profile considered here. The effect of the clad- 
ding is neglected, which is justified in view of the earlier 
given numerical values . Taking the average temperatures 
simplifies the mathematical solution without introducing 
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considerable errors. If more accuracy is desired, the obtained 
temperature values could always be corrected with a multipli- 
cation factor for point values or for optimal accuracy the 
temperature distribution model may be obtained by the solution 
of partial differential equations [3]. 



a single fuel element, whose diameter is small with respect to 
the core, may be written for the steady state heat transfer 
as : 



Under the above assumptions the heat balance for 




(III. 14) 




Fig. 7. Cross section of cylindrial 
fuel element. 
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The solution to this equation with the boundary condition 



dT 

dr 



r=0 



0 



T 



TF 



r=0 



is (see Appendix B) 



Q p = K A (TF - TAV) 



(III. 15) 



where K = Thermal conductivity, in this case an average 
value for the fuel cladding and coolant 



cal 



2 

sec cm °C 



2 

A = surface area of the fuel element [cm ] 



The thermal conductivity of the fuel and cladding and the 
physical properties of the coolant (density, viscosity, and 
specific heat) are considered to be constant and independent 
of temperature, which is a valid assumption for medium power 
variations of the reactor. For the above assumptions the 
product K . A may be considered constant, depending only upon 
the physical size of the fuel element. For the dynamic case 
of the heat transfer the heat capacitance of the fuel must be 
taken into account. Using the analogy of an electric circuit 
for heat flow, where heat flow is equivalent to current, 
temperature equivalent to voltage, and heat sink (dissipation) 
equivalent to a resistor, the following relationship can be 
written 




dTF 




(III. 16) 



dt 
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Coolant 






Fuel 




7T^77 




Fig. 8. Heat flow model. 



The definition of the heat capacitance is [7] 



A Q cal 

A T sec 



Therefore the above equation may be written as 

~ (TAV - TF) + A' N (III. 17) 

. 1 K A 

where — = 

T c f 

= fuel time constant 



cal 



cm 



sec cm 



cal 



sec 
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. _ A 



A' = 



= Proportionality factor 
b . Heat Removal 

The amount of heat removed by the coolant is given 



as 



Q = e * A, 



V 



AT 



(III. 18) 



where 3 = density of coolant 



gm 

3 

cm 



A - cross sectional area of coolant passage [citi ] 



V = coolant speed 



cm 

sec 



C = specific heat of coolant 



cal 



gm 



AT = coolant temperature rise within the element [°C] 
As shown earlier the temperature rise of the coolant is not 
a linear function of the position of a unit element within 
the coolant channel. A second nonlinearity accours in the 
coolant channel which is shown in Fig. 9. The temperature of 



Fuel 
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Fig. 9. Temperature in coolant path, 
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coolant at the walls of the fuel element are higher than at the 
center of the coolant channel. This nonlinearity may be 
neglected for the reasons stated earlier, and an average tem- 
perature profile will be assumed as shown in Fig. 9. The 
specific heat of the coolant and the coolant density will be 
assumed constant. The heat balance equation for the dynamic 
case is again derived with the aid of an electrical equivalent 
as shown in Fig. 10. 



Coolant 




Fig. 10. Electrical equivalent for heat removal 
by coolant from core . 



The heat flow equation is 



Q F 



C 



c 



dTAV 

dt 



+ Q 



AV 



(ix J. . ) 
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where Q ^ = Q (.from equation III. 18) 
which may also be written as 



Q 



AV 



where m 



m c AT 

m c (TCI - TH1) 
coolant flow rate 



gm 

sec 



After rearranging and using the relation 



TAV 



TCI + TH1 
2 



to eliminate TH1 , the fallowing equation is obtained 



dTAV 

dt 



1 TAV + ^ — TF + ^ — r TCI 
2 2 T 2 X 2 



(III. 20) 



where i— 

X 2 



KA r cal 2 "C 

tt ~ ~ cm — r 

C c [sec cm °C Ca 





feme gm °C cal ^ 1 

I C c sec cal gm °C sec 



2 . Heat Exchanger 

Assuming no heat loss to the outisde the following 
relation holds : The total change in heat of the secondary 

sodium within the heat exchanger. 

AQ X = AQ 2 (III. 21) 

where AQ^ = m c(TBIl - TB01) 

= the heat transferred from the primary sodium 
AQ 2 = KA (TBl - TB2) 

= the heat transferred to the secondary sodium 
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Using average temperature of primary and 
secondary side. 



Equation III. 21 becomes now 



KA (TBI - TB2 } = m c (TBI1 - TB01) 



(III. 22) 



or, after using the relation 



TBI 



TBI1+TB01 

2 



TBI1-TB2 = (TBIl-TBl) 



where Ml =. ^ P 



[ 



2 

gm cal sec cm °C 1 
sec gm °C cal 2 



cm 



-f dimensionless 



If now, as done before, the heat flow rate, heat capacitance 
of the primary sodium and the overall transfer coefficient is 
assumed constant, Ml may be seen as a heat exchanger constant. 

The dynamic equation for the secondary side is derived 
as in section II.B.l.b. 




where Q p = KA (TBI - TB2 ) 



Q s = m c (TH2 - TC2 ) 



C = heat capacitance of secondary sodium 

O 



or after rearranging 



dTB2 

dt 




TB2 + — TBI + 
T 3 



1 



TC2 



(III .23) 



T 



3 
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where 




1 = 2 m c f 1 

V c s L sec . 

3 . Steam Generator 

Since the flow of secondary sodium through the steam 
generator can be changed by the bypass valve in the secondary 
circuit the flow rate of sodium through the steam generator 
will be defined as x*m, where x is the valve opening 
(0 < x < 1) . 

x = 1 valve fully open (100% flow rate) 

x = 0 valve completely closed 
As in section II. B. 2. the following equation is obtained 



KA (TB3 - TS) = x m c (TBI2 - TB02) 



(III. 24) 



eliminating TB02 and rearranging gives 



TB3 - TS = — (TBI2 - TB3) 
M2 



where 



M2 



2 m c 
KA 



= a dimensionless steam generator constant 
under the assumptions of section II. B. 2. 
For the dynamic equation of the steamside, the principle used 
in section II.B.l.b. is again applied and the following 
equation obtained 



dT c 

q p " c s dF 



+ Q c 



(III. 25) 
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where Q p = KA CTB3 - TS) 

C g = heat capacitance of the steam 

Q s = heat consumed by the steam plant 
From the Carnot cycle it is known, that the work done in a 
thermodynamic machine is proportional to the difference of 
the added and rejected heat. 

W a <Q a - Q r ) 
but Q S = °A " °R 

and is proportional to the entrance temperature to the 
steam plant and Q D is proportional to the steam outlet tern- 
perature. Taking the average temperature it is seen that 

q ^ Ci TS Ci W 

from this relation Q g may be written as . 

Q c = K • \p • TS 

K = constant, assuming, constant flowrate, heat 

capacitance and linear temperature change with 

respect to pressure over the range of operation 

cal 

°C 

ip = dimensionless proportionality factor for the 
load demand (0 < ip < 1) 
ip = 1 ■+ 100% load 
\p = 0 -> 0 load 
equation (III. 25) is then 

KA(TB3 - TS) = Cg + K • ^ • TS 
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and after rearranging 



dTS 

dt 



where 




K' ip) TS + 




TB3 



dimensionless 



■] 



Cm. 26) 



C. HEAT TRANSPORT 

To obtain dynamic equations which include the time delays 
present within the heat transport system of the plant, the 
Taylor series expansion will be used. Throughout .this section 
it will be assumed, that there is no heat loss during the heat 
transport from one plant component to another. 

1 . Primary Circuit 

The transport delay, i.e., the time it takes for one 
unit volume element of sodium to get from the exit of the 
reactor (with temperature TH1) to the entrance of the heat 
exchanger (still with temperature TH1, but now called TBI1) is 
designated T4. Thus 
THl = TBI1 

but THl (t) = TBI1 (t + T 4 ) (III. 27) 

Expanding the right side of equation (III. 27) in a Taylor 
series 



TBIlCt+T 4 ) = TBIl(t) + T 4 dTB ^ . l t). + h.O.T. (III. 28) 



where H.O.T. = higher order terms. 
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Assuming slew temperature changes, which is reasonable 
for a large stationary power plant, all higher order terms may 
be neglected. 

Replacing TH1 and rearranging, equation (III. 27) may 
be written as 



dTBIl 

dt 



1 _ 

T 



TBI1 + 



T, 



TAV 




(III. 29) 



For the transport delay from the heat exchanger back to the 
reactor the same time duration will be assumed. 

Therefore TB01 = TCI 

but TB01(t) = TCI (t + T 4 ) (III. 30) 

Using again Taylor series expansion, replacing TBOl and 
rearranging, the dynamic equation for the heat transport from 
the heat exchanger back to the reactor is given as 



rlTf 1 ! 1 9 I 

= - ±- TCI + =r- TB1 - ±- TBI1 (III. 31) 

at 1 4 1 4 1 4 

2 . Secondary Loop 

The heat flow diagram for the secondary sodium loop is 
shown in Figure 11 where Q n (n= 1,2,. ..,5) is the rate at 
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Fig. 11. Heat flow diagram for secondary 
circuit. 



which heat is transported by a unit volume element. The 
function of the bypass valve X is to increase or decrease the 
amount of sodium flow through the steam generator: it does 
not influence the temperatures. This may be shown as follows: 

using the equivalent of an electric circuit, the "model 
equation" for node A is written: 

(III .32) 



Q 2 ~ Q 1 


" Q 3 




x m c 


TBI2 = m c 


TH2 - (1 - x)m c TH2 


x TBI2 


= TH2 - (1-x) 


TH2 


TBI2 


= TH2 
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Therefore the same procedure as in section III.C.l, may be 
used to arrive at a dynamic equation for the heat transport 
delay from the heat exchanger to the steam generator. The 
time duration of the transport delay will be called T^ . 

Since TH2 = TBI2 

and TH2 (t) = TBI2 (t + t 5 ) 

and, again using Taylor series expansion for TBI2(t + T,.) and 
replacing TH2 (t) gives 



dTBI2 (t) 
dt 



- TBI2 + =— TB2 
5 5 




(III. 33) 



For the heat transport from steam generator to heat exchange 
the "nodel equation" for "node" B is 



q 5 = q 3 



( TTT 
V J - 






4) 



m c TC2 = x m c TB02 + (1-x) m c TH2 

Assuming that the duration of the transport delay is again T^ 
for the flow through the bypass and for the flow from the 
steam generator, equation (III. 34) becomes 



TC2(t + T,) = x TB02 (t) + (1-x) TH2 (t) 

D 



(III. 35) 



Replacing TB02 and TH2 , using the relations 



M, 



x 



TB3 = 



TBI2 - TB3 
TB3 - TS 

TBI2 + TB02 



TB2 = 



TC2 + TII2 
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and expanding TC2 Ct + T^) in a Taylor series, equation (III. 35) 
becomes 



dTC2 

dt 



J~r TC2 + TB2 + ~~ TS + TBI2 (III. 36) 

1 5 5 5 5 



where 



T 



5 



T 



5 



2-x 



T 



n 

5 



T 



5 



2-2x 



T c (x+M2) . 

m III _ _2 

5 2x M2 



T c (x+M2) 

m iv _ _2 

i 5 x(x-M2) 
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IV. SYSTEM SIMULATION 



A. ANALOG COMPUTER PROGRAMMING 

The coefficients of the differential equations developed 
in section II were calculated for a steady state solution (see 
Appendix C) . Using the numerical values of these coefficients 
the following set of equations which describe the dynamic 
behaviour of the plant was obtained: 



D 


^ . 58 (N-D) 


TF 


= - .5 TF + .5 TAV + 725 N 


TAV 


= - .5333 TAV + .5 TF + 4.8333 TCI 


TCI 


== - .25 TCI - .0988 TBIl + .3488 TB2 


TBIl 


= - .25 TBIl + .5 TAV - .25 TCI 


TB2 


= -6.9626 TB2 + 3.0233 TBIl + 3.9394 TC2 


TBI2 


- - .25 TBI2 + .5 TB2 - .25 TC2 


TC2 


= - .2625 TC2 + .275 TS - .0375 TB2 + .025 TB2 


TS 


= - .0421 TS - .0168 • Tp • TS + .0421 TBI2 


N 


= D / (1-ICp) 



To simulate these equations on an analog computer and to obtain 
maximum accuracy, the gain coefficients must to be scaled with 
respect to the maximum voltage level of the analog computer 
(see Appendix D) . The scaled equation used for the simulation 
are given below. 



t 

D 


= .58 CN-D) 


TF 


= - .5 TF + .1833 TAV + .4833 N 


TAV 


= -5.3333 TAV + 1.3636 TV + 4.3939 TCI 
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TCI = - .25 TCI - .1235 TBI1 + .3314 TB2 
TBI1 = - .25 TBI1 + .44 TAV - .2 TCI 
TC2 = - .2625 TC2 + .2895 TS - .0474 TBI2 + .025 TB2 

TB2 = -6.9626 TB2 + 3. + 78 TBX1 + 3.9394 TC2 

TBI2 = - .25 TBI2 + .3958 TB2 - .1979 TC2 

TS = - .0421 TS - .0168 • 4 • TS + .0505 TBI2 

N = D/U-K^,) 

The temperature variables for example are no longer in °F but 
in volts . 

A block diagram, using both the symbols shown below and a 
linear flowgraph for representing the analog computer setup 
may now be drawn (Figure 13, 14). (See Appendix -E for the 
analog computer diagram corresponding to the flowgraphs shown 
in Figure 14) 



-t>- 

-o- 

o 




Adder Amplifier 
Integrator 
Potentiometer 
Multiplier 




Divider 



Figure 12 . Symbols for analog computer components 

It proved very valuable to write a simple Fortran executive 
program which performed the above calculations each time a 
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test or simulation run was made. This program contained the 
following steps: 

a) Read in steady state plant variable values, 
constants and chosen maximum plant variable 
values 

b) Calculate the numerical value for the coefficients 
of the differential equations 

c) Calculate magnitude scaling factors 

d) Calculate scaled coefficients for the analog 
computer components 

e) Calculate initial conditions for analog 
integrators 

f) Write out all calculated values and the scaled 
and unsealed differential equations 

g) Calculate the initial steady-state output values 
of the amplifiers 

h) Set all potentiometers 

i) Scan all amplifiers and compare with calculated 
values , write out errors 

The amount of time saved by using this executive program 
during the many runs performed was appreciable. The program 
had the following advantages : 

a) Equipment failures were detected immediately and 
isolated to only a few components by the test part of 
the program. 

b) The scaling procedure was done by the computer. Since 
the maximum plant variable values were not known, the 
optimum scaling coefficients had to be found by trial 
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and error in many runs . ^ By just changing one statement 
in the program (maximum variable value) the analog 
computer was ready for the next run within seconds. 

c) The documentation was done automatically by the digital 
computer. 

B. DYNAMIC CHECK 

A theoretical curve describing the behaviour of the reactor 
of the proposed plant following an accidental reactivity input 
was available from Ref. [2] . The validity of the dynamic 
equations describing the reactor was demonstrated with the aid 
of this curve. 

To check the solution of the analog computer the system 
equations were also solved on a digital computer (IBM 360/67) 
using a fourth order Runge Kutta routine with accuracy testing 
of single and double increment and double precision calculation. 
The results of these dynamic checks are stated below. 

1. Controlled Accidental Reactivity Input 

The accident considered assumes that two control rods 
drop into the core. This would correspond to a reactivity 
insertion of 100 $/sec. It is further assumed that scram 
Action takes place .2 sec. after the accident starts. 

A solution of the system equation to this accident was 
obtained on the digital computer with the above described 
routine . 

■^Because of the nonlinearities in the plant dynamic equa- 
tions, an analytical solution for the calculation of the 
maximum values cannot readily be obtained. 
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The digital computer solution gave exactly the same 
results for the temperature variables TAV and TF. The numerical 
values for the neutron density N (t) varied slightly in magni- 
tude and in the time of occurance. This slight difference can 
be explained by the simplifying assumptions which were made in 
deriving equation II. 5 and II. 6. The one group delayed 
neutron and the prompt jump approximation effects the time of 
occurance of the maximum value of N(t). By neglecting small 
reactivity coefficient (equation II. 10) the magnitude of N (t) 
is influenced. In spite of this small difference in the 
result, the equations used for the simulation of the reactor 
core are thought to be accurate, especially since, as stated 
above, the temperature TAV and TF followed very precisely the 
given data. 

2 . Change in Load 

The first part of the check on the accuracy of the 
analog computer solution was made by obtaining a digital compu- 
ter solution to a load change in the form of a step change of 
magnitude - 20% for 20 sec. The results are given in Table I 
and in Figure 15. 

3 . Change in Reactivity 

A second digital computer check of the analog computer 
solution consisted in introducing a disturbance in the reac- 
tivity input in the form of a +20 percent reactivity step 
change for 20 sec. The results are given in Table II and in 
Figure 16 . 
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Time (sec) TCI(.°F) TB2 C°F) TS (°F) 



0 


800; 2 


885.1 


750.2 


2 


801.1 


886.6 


762.3 


4 


803.9 


890.4 


773.1 


6 


808.2 


895.5 


782.9 


8 


813.3 


901.5 


791.9 


10 


820.2 


908.6 


800.2 


15 


835.3 


924.9 


818.5 


20 


840.1 


941.5 


834.7 


22 


844.6 


946.5 


827.6 


24 


852.2 


94 8.5 


821.7 


26 


857.4 


948.8 


816.9 


28 


858.3 


948.1 


812.9 


30 


860.6 


946.8 


809.2 


35 


855.2 


942.2 


804.0 


40 


850.4 


937.3 


800.0 



TABLE I 

Digital computer results ; temperature changes due to 
step change in load of -20% for 0 < t < 20 sec. 
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TS 






100 % 

80 % 



Fig. 15. Analog computer results Load change as in table I. 
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Time Csec) 


TCI C°F) 


TB2 C°F) 


TS C°F) 


0 


800.1 


885.1 


750.2 


2 


800.3 


887.8 


750.2 


4 


801.1 


891.6 


750.5 


6 


802.2 


893.9 


751.0 


8 


803.3 


895.5 


751.8 


10 


804.3 


896.8 


752.8 


15 


806.5 


899.4 


755.8 


20 


808.9 


901.8 


758.9 


22 


809.7 


900.6 


760.2 


24 


810.0 


897.9 


761.0 


26 


810.0 


896.4 


761.7 


28 


810.0 


895.6 


762.0 


30 


810.0 


895.3 


762.0 


35 


810.4 


895.2 


761.4 


40 


810 .6 


895.0 


760.5 



TABLE II 

Digital computer results; temperature changes due 
to step change in reactivity input of + 20 cents 
for 0 < t < 20 sec.- 
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TS 





760°F 

750°F 

TB2 

895 

885 




TCI 

810 

800 



React. 



3.36% 
3 . 5 6 % 



Fig. 16. Analog computer results Reactivity change as 
in table II. 
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4 . Discussion of the Dynamic Check Results 



It can be observed from the above results that the 
developed equations are valid for the plant under consideration. 
Comparison of the given digital computer solution and the 
analog computer graphs shows that the analog computer gives a 
solution which is within the bound of less then one percent 
exact (+ . 5%) . 

The results of the tests have demonstrated the validity 
of the simulation within the above bound of accuracy and under 
the earlier made assumptions . 
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VI . CONTROLLER 



A. CONTROL STRATEGY 

The inherent incapability of the plant to follow load 
changes which is due to the small reactivity feedback coef- 
ficient of the coolant, is a disadvantage in the control of a 
Fast Breeder Reactor. It makes complete control of the coolant 
temperature in the reactor necessary. An automatic controller 
for power level control must be developed whose two main 
purposes are 

a) coolant temperature level must be kept constant 

to obtain optimum breeding ratio and optimum perfor- 
mance of the reactor and 

b) steam conditions are to be kept constant for optimum 
performance of the steam plant. 

On the other hand several advantages which are inherent 
for the plant help the controlling action to achieve smooth 
and efficient power level changes. The reactor is automati- 
cally isolated from the load by the secondary sodium loop and 
to a certain degree by the primary sodium loop. Therefore, 
the inlet temperature to the core changes only very slowly. 
These two circuits also serve as a heat reservoir. This allows 
a temporarily larger rate of load change than the rate of power 
change which the reactor can achieve. 

The .rate of power change obtained from the reactor is 
limited by the maximum rate of change of reactivity. This 
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cent s 

value is taken as 2 /sec. 12] It is also limited by 

mechanical factors such as thermal stress within the reactor, 
heat exchanger or steam generator and speed of mechanically 
driven valves. The thermal stresses were not investigated 
here, it can, however, be stated that in none of the test runs 
the limiting values for the rate of change of temperature, 
e.g. the coolant in the core with 100°F/sec. , as given in 
Ref. [2] , was observed. The limiting factor in this investi- 
gation was found to be the rate of change of reactivity. 

The control strategy is based on the following thoughts: 

Because of the extreme long time constant, a pure pro- 
portional controller, which acts only following the occurance 
of errors in the to-be-controlled quantities was not sufficient. 
Therefore, an anticipatory or derivative control had to be 
implemented. To get the best derivative control action signal, 
the derivative or the rate of change of temperatures which had 
to lead to changes in the to-be-controlled temperature had to 
be sensed at the earliest possible time. This steady state - 
or at most very slow power level change - strategy is not 
valid for larger power level changes. In this case preprogrammed 
controlling action had to be implemented. 

B. CONTROLLER DESIGN 

The choice was given between a pure analog-logic controller 
or a digital controlling program. 

The analog controller had the advantage that it could have 
been implemented comparitively easy for a plant simulation of 
this size. On the other hand it was felt to be more realistic 
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to design an automatic controller which, uses the digital 
computer since modern power plants are equipped with digital 
computers used for many purposes - from bookkeeping of the 
plant datas to paycheck writing - that a digital controller 
would represent only a small and comparitively low cost 
addition to the main computer. Investigations have also shown 
[8] , that for large scale plant controllers digital controllers 
become more economical than analog controllers. Therefore, it 
was decided to design a digital controller. The following 
controlling devices were "built" in to the plant: 

a) control rods. It was assumed that the rate of change 
of reactivity could be varied. For the preprogrammed 
controlling action the maximum rate of 2 cents/sec. was 
assumed; for the steady state controlling action a 
value of . 1 times the maximum rate was found by 
experiment to give the smoothest control dynamics. 

b) Bypass valve in secondary circuit. The opening time 
(with maximum speed) for the preprogrammed phase was 
assumed to be 20 sec. from completely closed to fully 
open; for the small controlling action .2 of this value 
was used. 

Because of the constantly changing value of the bypass valve 
opening, several gain factors in the dynamic equations 
developed in section II changed during operation. As a result 
of this the analog computer setup had to be modified (see 
Appendix E) . 
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The above described controlling devices where controlled 
by the digital computer which in turn received the data from 
the simulated plant on the 




Fig. 17. Plant and control points. 



analog computer (see Fig. 17) . 

The control strategy called for 

a) preprogrammed power level change control. The control 
trajectory for this control was obtained by introducing 
load changes and correcting the errors in TCI and TS 
with manually operated potentiometer, which simulated 
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control rods and bypass valve opening. The resulting 
datas were evaluated with the aid of a least square 
straight line approximation program to give the 
control trajectory for the control reactivity and valve 
opening as a function of the load (see Appendix F) . 

b) Derivative control. 

Following the thought described in IV. 1, the rate of 

change in TAV and TB2 were calculated by storing values 
th 

every 10 of a second and then calculating the rate of 
change over one second. 

c) Proportional control. 

The to-be-controlled temperatures TCI and TS were 

th. 

compared every 10 of a second with the set value 
and the error was calculated. 

To take the addition of noise in actual plant measured 
variables into consideration, all measured values were put into 
a dead zone. This had also the advantage of smoothing the 
controlling action. 

A flow chart of the basic program for the digital controller 
is given in Figure 18. The individual parts of the program 
are described in detail in Appendix G. 
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Fig. 18. Flow chart for digital controller. 



53 



VI. RESULTS AND COMMENTS 



The results which were achieved by controlling the plant 
with the described controller are shown in this section as 
stripchart recordings in Figure 19 through Figure 21. The 
test inputs were step changes in load of different sizes and 
at different time intervals. 

Figure 19 shows the plant dynamics following a step change 
in load of -50%. This step change does not seem to be very 
realistic, it proves, however, that the initial overshoot in 
the to-be-controlled temperatures of little more than 100 °F 
is damped down very effectively to a secondary overshoot - or 
first undershoot - of about 20°F. The initial overshoot 
cannot be avoided since the controlling devices work at 
maximum allowable speed during the preprogrammed part. The 
same statements hold for the step change in load of +50% which 
was introduced after the transients had almost settled. Even 
a load change of this amount does not endanger the plant - due 
to the fast working and direct preprogrammed controller - which 
may be seen from the value of N which increases only about 12%. 

It is believed, that this test input and the resulting 
transient demonstrate the validity and necessity of the 
division of the control program into its two main parts, it 
also demonstrates effectiveness of the preprogrammed controller. 

Figure 20 shows the plant dynamics following a step 
change of different magnitudes and at time intervals which 
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are shorter than the settling time of the dynamics. A 20% 
step change in load produces only a 25°F, or 35°F respectively, 
overshoot in the controlled temperature. The increase in N 
is only 4%. The interplay of preprogrammed and small deviation 
control may also be seen from Figure 20 by considering the 
reactivity and bypass valve curve. 

Figure 21 was taken over a time of more than 1 hour. 

Though step changes in load of 10% do not seem to be very 
practical this recording shows a more realistic picture of the 
plant controller. It especially proves the accuracy of the 
steady state controller. After the settling of the dynamics 
this controller is able to keep the error in the controlled 
temperatures within a bound of less than + 3°F at a controlled 
temperature level of 800°F for TCI and 750°F for TS . 

The above results demonstrate the effectiveness of the 
designed controller for the plant layout simulated in this 
investigation. Since conventional stability calculations like 
root-locus plot, Bode plot or Nyquist diagram cannot be 
applied directly to a controller of the kind developed here, 
this topic is not included in this report. It can however be 
stated that no instabilities following any of the applied load 
changes were observed during the simulation. 

The addition of more plant components in the simulation 
of the system will, as stated in section II, not change the 
basic form of the dynamic system equation. It will however, 
change the dynamic response to a certain degree. The same 
is true for the nonlinearities which were not taken into account 
in the simulation in this report. 
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Further investigations can include the above considera- 
tions, however, a simulation of a complete plant which includes 
all components and takes all nonlinearities into account is 
only possible on a very large scale digital, or better, hybrid 
computer. If a computer of this size is not available an 
accurate simulation can be done component by component as in 
Ref. [3] for a LMFBR. Accurate stability investigation may 
then be performed by simulation. 
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10 min 0 Time (. 1 min/sec.) 

Fig. 19. Loadchange (50% Step) with control action. 
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10 min 0 Time (.1 min. /sec.) 
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20 min 10 min 0 Time (. 1 min. /sec.) 

Fig. 20. Load changes (Small steps) with control action. 
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Fig. 21. Longterm control action. 
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APPENDIX A 



A. CALCULATION OF THE REACTIVITY FEEDBACK COEFFICIENTS 
The reactivity coefficients considered are 

a) Doppler coefficient 

b) Fuel expansion coefficient 

a) The Doppler coefficient is a nonlinear function 
of the fuel temperature and is shown in Fig. 22. 




Fig. 22. Doppler coefficient. 



Since the operating temperature of the reactor is at 2400 °F an 

'-5 

average doppler coefficient of -.35*10 could be taken. It 
can be seen, that even at temperature variations of +400°, 
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which, have to be avoided in a fast reactor, the change in 

-5 

doppler coefficient is at most +_.5*10 . To account for a 

certain degrees of uncertainty a conservative numerical value 
-5 

of -.3-10 will be used. Changing this value to dollar units 
. . -4 

of reactivity with (3 taken as 29.88x10 the doppler coefficient 
is obtained as 



K d = — i- . 3xlQ = -1 . 04xl0’~ 3 dollars 

29 . 88xl0~ 4 

b) The fuel expansion coefficient is considerably 

smaller than the doppler coefficient and since its reliability 

for axial fuels is highly uncertain the conservative numerical 

—6 

value of -1.0*10 will be taken. The conversion into dollar 
units reactivity gives: 



k exp 



- 1.-10 6 
29 . 88x10" 



.34x10 ^ dollars 



The total feedback reactivity is therefore given as 

k f = k d + k ex 

= -.00104 - .00034 
= -.00138 
= -.0014 dollars 
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APPENDIX B 



A. SOLUTION OF THE STEADY STATE, SECOND ORDER DIFFERENTIAL 
EQUATION FOR THE RADIAL HEAT FLOW IN A CYLINDRICAL FUEL 
ELEMENT 

The heat balance equation is given as 



Q 2 nr dr L = Q . , - Q 

r+dr r 



The differential heat quantities 



Q = Q + dQ = q + -5— ; = 
r+dr r r r dr 



dr 



(B.l) 



(b . 2 ) 



dT 

Q = - K, A 

r f dr 



T , ~ T dT 
- K, 2irr L — 
f dr 



(B . 3) 



equation (B.2) becomes now 

? 

OTV dT d“T dT- , 

r+dr r dr ,2 dr 

dr 



(B . 4 ) 



or 



Qp 2nr dr L 



- 2 it L K(r + ~) dr 

dr 2 dr 



(B . 5) 



which, after concellation and rearranging, may be written as 



K a-| + £ ® + 0 

2 r dr F 
dr 



(B . 6 ) 



This differential equation can be solved by standard methods, 

the solution is 
2 



T= - Q F4K +C 1 £nr+C 2 



(B . 7) 
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The boundary conditions are 



dT 

dr 



r =0 



0 



T 



r =0 




T I = T 
' r=R c 



Applying these boundary conditions gives : 



T - T 
c s 




Since 



Q R = nr L Q f 



(B. 8 ) 



it follow that 



Qg — ttR L Qp 



(B . 9 ) 



or, combining (B. 8 ) and (B.9) 



Q s = 4 tt K L(T c - T s ) 



(B.10) 



Setting 2u R L = A, the surface of the fuel element 



Q s = 2 K f A 



T C " T S 
R 



(B.ll) 



Equation (B.ll) holds for the region inside the fuel. For 
the cladding region K ^ 3 would have to be used instead of 
and for the region in the coolant K ^ an ^. Since average 
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temperature are used and in addition the individual thermal 
conductivities are considered constant, an average thermal 
conductivity will be taken and (B.3) may be written as 

Q = - K 2-rrr L (B.12) 

r dr 

where r = the distance measured from the center of the fuel 
element. 

T c is taken as TF, the average fuel temperature and T g 
as TAV, the average coolant temperature. 

Equation (B.5) beconjes then 

Q„ = KA (TF - TAV) (B.13) 

r 
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APPENDIX C 



A. CALCULATION OF THE GAIN COEFFICIENTS OF THE SYSTEMS 
DIFFERENTIAL EQUATIONS 

The gain coefficients in the individual equations may be 
solved for by using the numerical values given for a specific 
plant. This would require a certain amount of engineering 
approximation, since exact values are not available for all 
quantities. Instead of following this approach, only given 
steady state temperature levels for 100% power (\p =1) at a 
bypass valve opening of 95% (X = .95) and given time constants 
were used. These values are {9] : 
a) Steady state temperatures 



TF 


= 


2400 °F 


TH1 


- 


1100°F = TBI1 


TCI 


= 


800 °F = TB01 


Th2 


“ 


1050°F = TBI2 


TC2 


- 


720 °F = TB02 


TS 


= 


750 °F 


b) Time 


constants 


T1 


= 


2 sec = T2 


T3 


= 


4 sec = T4 = 


T5 


- 


10 sec 


T7 


= 


. 1 sec 


Additional constants used [3] 


X 


= 


,58 


B 


— 


.002988 
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Setting these constants into the obtained dynamic equation 

these equations were "forced" into an accurate steady state 

solution by solving for the unknown gain coefficient in the 

dT 

steady state form of the differential equation (~ = 0) . 

This procedure is demonstrated with one example: 

The equation for heat transfer within the core was obtained 
as : 



TF 



TF + TAV + A'N 
1 1 



which becomes in its steady state form 
0 = - | (2400) + j (950) + A' (1) 



and solving for A s gives 



A' = 750 

which is the desired coefficient for N. Using this approach 
for all equations individually leads to the system equations 
given in section II. 
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APPENDIX D 



A. SCALING OF THE SYSTEM EQUATIONS 

As mentioned in section III the actual equations must be 
scaled with respect to the maximum voltage levels of the 
analog computer. It is desirable to get maximum output values 
of the equation variables close to the computer maximum 
voltage level to maintain maximum accuracy in the solution of 
the differential equations. To be able to perform this scaling 
in an optimum way an exact knowledge of the maximum values 
of the plant variables is necessary. If these variables are 
not known exactly, optimum scaling is done by trial and error: 
a maximum value is assumed, the equations are scaled for a 
computer run and the output signal is observed. If this 
output signal is either to high or too low, the equations have 
to be rescaled. 

The magnitude scaling procedure uses scaling coefficients 
which convert given plant variables with their individual 
dimensions to computer variables with the deminsion volt. 

Taking the example of TAV it is written 



TAV 

where SLTAV = 



TAV = 



SLTAV • TAV 



scaling coefficient for the plant variable 



TAV with dimension 



V 



computer variable for average coolant temper- 
ature in core with dimension (V) 
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The analog computer used here (COMCOR Ci 5000) had a 
maximum voltage level of + 100V. Assuming a maximum value of 
1100°F for TAV the scaling coefficient is obtained as 



SLTAV = 


TAVX _ 1100 . . [°F 

TAVX 100 11 [v _ 



Using this procedure the following numerical values for the 
individual scaling coefficient proved to be best suited: 



SLD 


= .02 


SLTF 


= 30. 


SLTAV 


= 11. 


SLTCl 


I! 

H 1 

O 

• 


STBI1 


= 12.5 


SLTB2 


= 9.5 


STBI2 


= 12. 


SLTC2 


= 9.5 


SLTS 


= 10. 



The scaling procedure for a single dynamic equation 
proceeds as follows : 

Taking equation (11.20) 

TAV = 5.3333 TAV + .5 TF + 4.8333 TCI 

Replacing the plant variables by computer variables gives 



SLTAV • 


TAV = - 5.3333* SLTAV • TAV + .5*SLTF • TF 




+ 4.833 - SLTCl • TCI 



Dividing through by SLTAV and using numerical values gives 
TAV = - 5.3333 TAV + 1.3636 TF + 4.3939 TCl 
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Using this procedure for all equations and steady state 
values gives the set of scaled equation given in III and the 
corresponding initial conditions for the integrators of the 
analog computer. These initial conditions were found to be 



TF 


(0) 


= 


• 

O 

00 


(V) 


TAV 


(0) 


= 


86.3636 


(V) 


TCL 


(0) 


= 


80 


Cv) 


TBI1 


(0) 


= 


88 


(V) 


TB2 


(0) 


= 


93.1579 


CV) 


TBI2 


(0) 


= 


87.5 


(V) 


TC2 


(0) 


= 


75.7895 


(V) 


TS 


(0) 


= 


75. 


CV) 



The steady state control rod reactivity required is found 
to be 3.36 dollar which, after being scaled, corresponds to 
33.6 (V) . 
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APPENDIX E 



A. ANALOG COMPUTER SETUP FOR CONTROLLER IMPLEMENTATION 

It was shown in II, that several gain factors in the system 
dynamic equation were a function of the bypass valve opening 
X. Since the bypass is used as one of the controlling ele- 
ments in the plant and therefore continually changes its 
numerical value, a way had to be found to continually change 
the gain factors which are functions of X. For the simulation 
without controller the gain factors were implemented as com- 
puter set potentiometers, which take a considerable amount of 
.time to be set or adjusted. A continuous dynamic solution or 
simulation of the plant was therefore not possible with this 
setup . 

The following method was used to get continuously changing 
gain factors and therefore a continuous simulation: 

The potentiometers were replaced by electronic multipliers 
which had as input the to be multiplied variable (temperature) 
and the output of a trunkline which was controlled by the 
digital computer. Figure 23 




Fig. 23. Implementation of varying gain 
factors . 
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shows the analog implementation for the second gain factor 
in the differential equation for TS (K**JwTS) . The voltage 
representing TS is multiplied by a voltage representing rp 
(0 [V] < ip < 100 IV] corresponds to 0 < ip < 100% load) . ip 
was changed by a manually set potentiometer on the analog 
board. The output of this multiplier is fed into another 
multiplier with the trunkline as the second input. Each 
time (scan intervall = .1 sec) the controller changed X, it 
simultaneously calculated all dependent gain factors and 
adjusted the corresponding trunkline. Due to the electronic 
multiplication continuous adjustment of the gain factors with- 
out any timelag was thereby achieved. 

Figure 24 shows the resulting analog computer setup for 
the plant part without the reactor. For more clarity of the 
graph no numerical values were set in and all initial condi- 
tion inputs to the integrators with the corresponding 
potentiometers were left out. Figure 24 corresponds to the 
linear flowgraph of the plant system equation given in III. 

Since there are no changing gain factors in the reactor 
part of the simulation Figure 13 in III is not changed. 
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tf 



-TAV 
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APPENDIX F 



A. CONTROL TRAJECTORY FOR PREPROGRAMMED POWER LEVEL 

CHANGE CONTROL 

300 data points were obtained by changing the load over 
the whole range from 100% to 0% in small steps. For these 
runs the analog computer was in 1000 : 1 timescale mode, 
therefore the steady state error for a certain load was 
immediately available. (see Figure 27). Using hand set 
potentiometers voltages 'were supplied to 

a) the analog computer which corresponded to control 
reactivity 

b) the digital computer which corresponded to the 
bypass valve position. 

The digital computer calculated the variable gain factors and 
adjusted the trunlines. 

The "control reactivity voltage" and the "bypass voltage" 
were adjusted until the steady state error was zero. The 
obtained data points were then used to obtain, with the aid 
of a least square polonomial fit program, straight line 
approximations as precomputed control trajectory for load 
changes in the plant. These data points are shown in Figure 
25 and Figure 26 . 

As may be seen from Figure 27, it is, with the controlling 
aids of this plant simulation possible to control TCI 
completely over the whole load range, TS can only be controlled 
over a load range from 100% load to about 20% load. The steady 
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state error for TS at zero load is about +80°F. This is not 
a restriction on the controllability of the plant. 
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Fig. 25. Experimental data points . 

Control Trajectory of reactivity 
Scale X: 1 inch = 10% load 

Y: 1 inch = dollar reactivity 
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Fig. 26. Experimental Data points. 

Control Trajectory of bypass opening 
Scale: X: 1 inch = 10% load 

Y: 1 inch = 10% valve opening 
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APPENDIX G 



A. DIGITAL CONTROL PROGRAM 

1. General 

The scan interval was chosen as .1 sec., i.e., the 
digital computer went through the designed control loop every 
.1 sec. and performed the necessary control action by scanning 
the assigned variables, performing the necessary calculations 
to get the derivatives or errors, and then going through all 
the logic steps, implemented as Fortran IF-statements . If 
control action was necessary it was performed by calling the 
following subroutines. 

a) REACP : Increase reactivity at the rate of 2 cents/ 

sec. Because of the scan interval the voltage of the 
trunkline which represented reactivity input was 
increased by an amount which corresponded to .2 cents/ 
.1 sec. Inspite of being a step change the short 
time interval in which these changes were performed 
made the controlling action an almost continuous one. 
The time delay of an electric motor which would be 
used to perform this controlling action in a real 
plant, was neglected. This is justifiable in view 

of the short time constant of a motor with respect 
to the time constants of the plant. 

b) REACM: Decrease reactivity at a rate of 2 cents/sec. 
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APPENDIX G 



A. DIGITAL CONTROL PROGRAM 

1 . General 

The scan interval was chosen as .1 sec., i.e., the 
digital computer went through the designed control loop every 
.1 sec. and performed the necessary control action by scanning 
the assigned variables, performing the necessary calculations 
to get the derivatives or errors, and then going through all 
the logic steps, implemented as Fortran IF-statements . If 
control action was necessary it was performed by calling the 
following subroutines. 

a) REACP : Increase reactivity at the rate of 2 cents/ 

sec. Because of the scan interval the voltage of the 
trunkline which represented reactivity input was 
increased by an amount which corresponded to .2 cents/ 
.1 sec. Inspite of being a step change the short 
time interval in which these changes were performed 
made the controlling action an almost continuous one. 
The time delay of an electric motor which would be 
used to perform this controlling action in a real 
plant, was neglected. This is justifiable in view 

of the short time constant of a motor with respect 
to the time constants of the plant. 

b) REACM: Decrease reactivity at a rate of 2 cents/sec. 
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c) RACP and RACM: same as under a) and b) but with 

the rate divided by 10. 

Subroutine REACP and REACM were called if controlling 
action was performed by preprogrammed controller. 
Subroutine RACP and RACM were called if controlling 
action was performed by small deviation controller. 

d) VENTP: Open bypass valve at described rate (0 - 100% 

in 20 sec.) . The subroutine changed the value of X, 
calculated the new T gain factor and adjusted the 
corresponding trunklines to the electronic multipliers. 
The same statements made in a) with respect to 
continuous action and time constant hold also here. 

e) VENTM: Decrease bypass valve opening at the described 

rate . 

f) VETP and VETM: same as under d) and e) but with 

rate divided by 5. 

Subroutine VENTP and VENTM were called by the pre- 
programmed controller subroutine VETP and VETM by the 
small deviation controller. 

The digital controller program was divided in 3 parts : 

a) Routing 

b) Preprogrammed controller 

c) Small deviation or steady state controller. 

These 3 parts are described below. 

2 . Routing 

The decision whether to go into the preprogrammed 
part of the controller or to the small deviation controller 
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was based on the rate of change of the load. For a load 
change of less than . 2%/sec. the control action was performed 
by the small deviation controller. The flowchart of the 
routing part is given in Fig. 28. 

3 . Preprogrammed controller 

The purpose of this controller was, to bring the 
controlling devices directly to their set points, i.e. to the 
points which were calculated by the straight line approximation 
- as a function of the load - to lie on the control trajectory. 
This controller performed its function without regard of the 
dynamics of the plant. The load was scanned, the set point 
calculated and the respective subprogram (Controlling action) 
called. Fig. 29 shows a flow chart of this part of the 
digital controller. 

4 . Small Deviation or Steady State Controller 

As stated in III this part of the controller used 
proportional and derivative control. The decisions of the 
digital computer were based on the logic tables given as 
Table III and Table IV. The procedure used here is the same 
as in analog logic circuit, instead of writing expressions in 
Bolean algebra, Fortran IF-statements were written for the 
logic decisions of the digital computer. 

Controlling action for change of reactivity was based 
on TAV and ETC1 (error in TCI) . The logic table used is 
Table III. 

Controlling action for change of bypass valve value 
was based on TB2 and ETS . The logic table used is Table IV. 
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With, respect to earlier thoughts any decision made on 
the basis of Table III was weighted higher and if controlling 
action was required, Table IV was neglected, i.e. the controller 
tried first to bring ETCI to zero and keep it zero before ETS 
was controlled. This was done because of the dynamic of the 
plant. No logic scheme could be arrived at, which would have 
made combined control possible. It will be shown in section 
V, that the performance of this controller was very good. A 
flow diagram of this part of the digital controller is given 
in Figure 30 . 
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Fig. 28. Routing part of digital controller. 
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Fig. 29. Preprogrammed controller part of digital controller. 
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Logic Tables for Small Deviation Controller 



TAV 

0 

+ 

0 



ETC1 

0 

0 

0 

+ 



+ + 
+ 

0 

+ 



Reactiv . 

•/. 

D 

I 

D 

D 

•/. 

I 

•/. 

I 



Table III. Reactivity control. 



TB2 

0 

+ 

0 

+ 

0 

+ 



ETS 

0 

0 

0 

+ 

+ 

+ 



B. P. valve 

• /• 

D 

I 

D 

D 

I 

I 



Table IV. Bypass valve control. 



Symbols : 

+ Error or Deriv. > 0 

Error or Deriv. < 0 

0 Error or Driv. = 0 

D Decrease 

1 Increase 

./. No action 
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Fig. 30. Small deviation controller part of 
digital controller. 
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APPENDIX H 



A. DIGITAL COMPUTER PROGRAM 

The following pages are a copy of the digital computer 
program which was used for the complete Hybrid computer 
simulation and controller implementation. The program includes 
all parts described in the earlier sections and appendixes. 
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CH9SSEN VARIABLES 



calculated initial conditions 
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tsx=sx 
T $ ^ = • 1 * S P 



TS'3 = SG 
A K ” = F V A L < 38 ) 

CALL DAC (l/TSX/?/TSS/3j»TSR/4/TSYj>biTSA'i6/TSG/8/TSP) 
CALL DAC(9/ A<R/10/X) 
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AL q AO( I ) = ALB AD ( I + l ) 
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